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A new single lane car following model of traffic flow is pre- 
sented. The model is inertial and free of collisions, ft demon- 
strates experimentally observed features of traffic flow such 
as the existence of three regimes: free, fiuctuative (synchro- 
nized) and congested (jammed) flow; bistability of free and 
fiuctuative states in a certain range of densities, which causes 
the hysteresis in transitions between these states; jumps in 
the density-flux plane in the fiuctuative regime and gradual 
spatial transition from synchronized to free flow. Our model 
suggests that in the fiuctuative regime there exist many stable 
states with different wavelengths, and that the velocity fluc- 
tuations in the congested flow regime decay approximately 
according to a power law in time. 



In the last years, growing effort has been made in un- 
derstanding traffic flow dynamics. Recent experiments 
show that traffic flow demonstrates complex physi- 
cal phenomena, among which are: 

(i) The existence of three states: free flow, "synchro- 
nized" (or "fiuctuative") flow and traffic jams (for 
low, intermediate and high densities correspond- 
ingly). The second state has two essential features: 
synchronization of flow in different lanes (for the 
multilane traffic) and fluctuation performed by the 
system in density-flux plane. Since our model is 
single-lane we will refer to this state as "fiuctua- 
tive" . 

(ii) Hysteresis which is observed in transitions between 
the free and the fiuctuative flow. 

(iii) Long survival time of traffic jams. 

Modeling of traffic fiow is traditionally performed using 
two approaches. The microscopic, or car-following mod- 
els approach, which describes the nearest-neighbor inter- 
action between two consecutive cars and investigates its 
infiuence on the fiow (see e.g. [3-5]), and the macro- 
scopic, or continuous models approach, which represents 
the fiowing traffic as a continuous media and describes it 
using the hydrodynamical partial differential equations 
(see e.g. [7-9]). Wide surveys of these models are given 
in [10-12]. 

In this Letter we inroduce an inertial single lane car- 
following model, which is free of collisions. We study the 



model both numerically and analytically and find the ex- 
istence of three regimes in traffic flow: free flow regime at 
low densities (where each car moves with almost a con- 
stant velocity), fiuctuative fiow regime at inetrmediate 
densities (where stable periodic oscillations of velocities 
of all cars are observed) and congested or jammed fiow 
regime at high densities (where due to high density all 
the cars tend to move with the same, relatively small 
velocity) . Our model predicts the existence of many in- 
homogeneous stable states in the fiuctuative regime and 
demonstrates hysteresis in transitions between free and 
fiuctuative regimes. The experimentally observed long 
survival time of jams may be explained by our finding 
that the fiuctuations in the congested fiow regime decay 
slowly according to a power law. 

To formulate the model we assume that car accelera- 
tion is affected by three factors: 

(a) aspiration to keep safety time gap from the car 
ahead, 

(b) pre-braking if the car ahead is much slower, 

(c) aspiration not to exceed significantly the permitted 
velocity. 

In mathematical description, the acceleration of the nth 
car a„ is given by a sum of three terms depending on 
its coordinate a;„, velocity v„, distance to the car ahead 
Axn = Xn+i — Xn and the velocities difference Aw„ = 
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where ^ is a sensitivity parameter, D is the minimal 
distance between consecutive cars, Vper is the permitted 
velocity and A: is a constant. The safety distance Aa:° = 
VnT + D depends on T, which is the safety time gap 
constant. The function Z is defined as Z{x) = (x+|a;|)/2. 
Note that Eq.(§) can be generalized by adding a noise 
term. 

In the following we discuss in more details the terms 
in the right side of (|l|): 

(a) The first term plays an important role when ve- 
locity difference between consequtive cars is rela- 
tively small. In this case the nth car accelerates if 
Axn > Ax'^ and brakes if Axn < Ax'^. 

The choice of function in this term is not unique. 
Replacing it by other functions of Axn which are 
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increasing, equal to zero if Aa;„ = and tend 
to —oo if Axn 0, such as Alog(Aa;„/Aa;°), gives 
similar results. 

(b) The second term plays an important role when 
Vn ^ Vn+1- A car getting close to a much slower 
car starts braking even if Axn > Ax^. This term 
corresponds to the negative acceleration needed to 
reduce |Ai;„| to as Ax„ D. 

(c) The dissipative third term is a repulsive force acting 
when the velocity exceeds the permitted velocity. 

Unlike optimal velocity models the acceleration in 
our model depends explicitly on Ax which enables us to 
make the flow free of collisions. 

The motion of cars is therefore described by the fol- 
lowing system of ordinary differential equations 



A{1 



^n + l—Xji ' 



(2) 



kZ(v„ 



n — 1,2, ... N with periodic boundary conditions 

N 

xn+i=xi-\ , Vn+1 = Vi. 

p 

A solution of Eqs. (||) which corresponds to homoge- 
neous flow is 



A(l~Dp)+kVp^ 
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In the following numerical results we use parameters 
values Vp^r = 25{m/s), T = 2(s), D ^ 5(to), 1 < ^ < 
5(m/s2) and k = 2(s-i). 

The flux-density relation (often called the fundamental 
diagram) for the homogeneous flow is shown in Fig. 1(a) 
as a dashed line. Comparison of this curve with the 
fundamental diagrams (solid lines) obtained by the nu- 
merical solution of equations (j^) for different values of 
A starting from nonhomogeneous initial conditions indi- 
cates that for values of p smaller than some critical value 
pi or greater than another critical value p2 the flux is 
the same, while for the intermediate values of density 
{pi < p < P2) the measured flux is considerably lower 
than the homogeneous solution flux. Plotting the vari- 

N 

ance of velocities = ("^n ~ ("f^))^]^^^ (where (u) 

ri=l 

is the average velocity) against p (Fig. 1(b)) shows the 



existence of velocity fluctuations for pi < p < p2. We 
can therefore define three regimes in traffic flow: the 
free flow regime (p < pi), the fluctuative flow regime 
{pi < p < P2) and the congested flow regime {p > p2). 
Note that the flow in the first and the last regimes is 
homogeneous. Note also that for small values of A p2 is 
greater than the maximal possible density Pmax = ^/D 
and the congested fiow regime does not exist. See Fig. 
1(b) for A = 2. This finding is supported by the analyt- 
ical results shown below. 




FIG. 1. (a) Fundamental diagram for A = 5,3, 2(m/s^) 
(top to bottom). Dashed line corresponds to the homoge- 
neous solution, (b) Ratio of variance of velocities to the av- 
erage velocity for A = 2, 3, 5(m/s^) (top to bottom), (c) 
Qualitative plot of function S[p). (d),(e) Hysteresis loops in 
transitions between free and fluctuative flow states for A = 3 
, arrows show the direction of changing the global density, 
(f) Results of local measurements of density and flux in free 
(almost straight line) and fluctuative regimes. 

In order to estimate the values of pi and p2 we analyse 
the stability of the homogeneous flow solution. The lin- 
earization of Eqs. (^) near the homogeneuos flow solution 
(^) in variables ^„ — Xn ^ x'^ has the form 



(4) 



AT+kTv 



where p = ATp + k, q = ATp+k — 

^ and p — ATp, q — Ap otherwise 



kD 



Ap^ for p < 



D+Tv. 

A solution of equation 



can be written as 
= exp{jan -|- zt}. 



(5) 
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where a = (k = 0, . . . , N — 1) and z - a complex 
number. Substituting into (||) we obtain the algebraic 
equation for z 



+ pz-q{e"^ - 1) = 0. 



(6) 



Each of the N equations (||) has two solutions. These 
2N different complex numbers are the eigenvalues of sys- 
tem (|4|). One of them (which corresponds to k = 0) is 
equal to zero regardless of values of parameters. In this 
case all ^„ in (||) are equal to a constant and belong to 
the one-dimensional subspace of equilibria of system (Q) 
(defined by equations £.1 — ■ ■ ■ = , — ■ ■ ■ — £,n — 0) ■ 
This indicates that the disturbed state Xn for z = is 
also homogeneous. For z 7^ ^„ in (||) is a wave with 
increasing or decreasing amplitude. Therefore, if we find 
conditions under which other 2N — 1 eigenvalues have 
negative real parts (the magnitude of wave (|^) decreases 
with time) we can say that under these conditions the 
homogeneous flow solution (0) is stable. 

Following the approach oflq] we can derive this condi- 
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tion as Y > 2 or S{p) > 2, where 

{ {ATp+kf ^ 1 

p-'A(AT+kv„„T+kD) ' I'' — D+Tv„^^' 
(7) 

A qualitative plot of S{p) is sketched in Fig. 1(c). From 
this figure it follows that depending on p we have three 
regimes of stability/instability of the homogeneous flow 
solution, li p < p' (free flow) or p > p" (congested flow) 
the homogeneous flow solution is stable and \i p' < p < p" 
it is unstable, where p' — jqjf^ — and p" = Note 
that there are possible sets of parameters under which 
the minimum of the left part of S{p) can be less than 
2 and the flow can have flve different regimes of stabil- 
ity/instability. Nevertheless, under the set of parameters 
specified above we have up to three regimes, where the 
third regime does not exist for p" > pmax [A < 2D /T"^). 



Our numerical simulations show that p2 



but p' 



is considerably greater than pi, thus we expect that for 
pi < p < p' both homogeneous and fluctuative states are 
stable. 

In the fluctuative regime (pi < p < P2) the flow is 
characterized by presence of humps (dense regions) mov- 
ing backwards or forwards. When the flow has stabilized 
the humps are equidistant and the evolution of traffic in 
time and space resembles the spreading of a wave. The 
existence of a fiuctuative regime was predicted by other 
car-following (e.g. ||^ where it was called "jammed fiow") 
and continuous (e.g. |9j, where it was called "recurring 
humps state" ) models and measured experimentally [Q . 

Simulations of our model show that the fiuctuative fiow 
state is not unique. Figs. 2(a-c) present the cars veloc- 
ities after the fiuctuative fiow regime has stabilized for 
three different initial conditions. It can be seen that the 



"wavelengths" of these states are different. Fig. 2(d) 
presents the convergence of fiux in these experiments to 
distinct values. Our simulations also show the existence 
of solutions with other "wavelengths" and flux values. 
Fig. 2(e) shows the fundamental diargams for three dif- 
ferent wavelengths. Consequently, depending on initial 
conditions different stable fiuctuative states emerge with 
different values of flux and distances between neighbor- 
ing humps. This indicates that for pi < p < P2 system 
(^) has many stable periodic (in Axn , Vn variables) solu- 
tions, and hence in the 27V-dimensional space of variables 
Aa;„, Vn there exist many attractive limit cycles. 
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FIG. 2. Three different stable states in the fluctuative 
regime, obtained from different initial conditions. Global den- 
sity p = 0.06(«e/i/m), yl = 3(m/s^). (a-c) Cars velocities, (d) 
Convergence of flux to different values in these three exper- 
iments, (e) Fundamental diagrams for three different stable 
fluctuative states with wavelengths 20, 5 and 6.67 cars (top 
to bottom). A dashed line corresponds to the homogeneous 
solution. 

As follows from above for pi < p < p' not only fiuctu- 
ative flow solutions are stable, but also the homogeneous 
flow solution. This bistability is the origin of hysteresis 
in transitions between free and fiuctuative fiow regimes. 
Such bistability was observed experimentally [Q and was 
found in other models [||j^J^. Fig. 1(d) shows a hys- 
teresis loop in the density-flux plane. The upper curve is 
obtained by increasing the density of cars adiabatically 
preserving the road length L [Q . It can be seen that up 
to the value of density p' the homogeneous flow is pre- 
served. The lower curve was obtained by adiabatically 
decreasing the density in the same manner. While de- 
creasing the density the flow remains fluctuative even for 
p < p' . Fig. 1(e) presents the hysteresis loop in the global 
density - velocities fluctuations plane. 

Our results also illustrate the well-known phenomenon 
||l|,||j|J|,|l^ of jumps which the system performs in the 
density-flux plane in the fluctuative flow regime when 
the density and the flux are measured locally. In our 
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numerical simulation (Fig. 1(f)) we started from a value 
of density below p', increased it gradually in the described 
above manner up to a value greater than p' and decreased 
it back. These jumps may be explained by our finding of 
many stable states in the fluctuative regime. 

Our model also demonstrates the gradual spatial tran- 
sition from the fluctuative to free flow in the downstream 
direction which was measured by . The results of local 
measurements of density and flux at different distances 
from an on-rump p5[ | are shown in Fig. 3. which is in 
good agreement with Fig. 3 of j^]. 



FIG. 4. Decay of the velocities variance in time in the con- 
gested flow regime. Road length is 1, 4 and 16 km (bottom 
to top). The dashed line has a slope -0.21 

In summary, we present a single lane car-following 
model which explains important features of traffic ob- 
served experimentally. The model predicts the existence 
of many stable periodic states in the fluctuative (synchro- 
nized) flow regime. 

We wish to thank S. Schwarzer for usefull discussion. 
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FIG. 3. Results of local measurements of flux and density 
at different distances from the on-ramp. (a) 500 m upstream 
the on-rump, (b) and (c) 250 and 1000 m downstream respec- 
tively 

In the congested flow regime the only stable solution 
is the homogeneous flow solution. We have not found 
evidence of existence of bistability or hysteresis in transi- 
tions between the fluctuative and congested flow regimes. 
Starting from random initial conditions, we observe that 
initial fluctuations of the velocity seem to decay accord- 
ing to a power law 



p-t/r 



t > t* 



(8) 



where t* is the crossover time between the power law 
and exponential decay. We find t* ~ and t ^ 
with z = 2.0 ±0.1. These results are qualitivly similar to 
that obtained by |l|] for a cellular automata model 
but with different values of exponents. The result z k, 2 
seems to be in agreement with random walk arguments 
of For the parameters values A = A, p = 0.15 we 
get /3 « 0.21 ± 0.04 (Fig.4). 
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